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We demonstrate the new specific phenomenon of the long-time resonant energy exchange in the 
carbon nanotubes (CNTs) with the radial breathing mode as a carrier. This process turns out 
to be stable in the wide range of the excitation energy. It is shown that the modified nonlinear 
Schrodinger equation, obtained in the framework of nonlinear elastic thin shell theory, describes the 
CNT nonlinear dynamics in the considered frequency band. The numerical integration of the thin 
shell theory equations confirms the results of the analytical study. 
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I. INTRODUCTION 

The study of nonlinear dynamics of the carbon nan¬ 
otubes (CNTs) is neccessary to describe the large elastic 
deformation of the CNTsi^, the processes of the phonon- 
phonon interactions^ — and the thermal properties of the 
nanotubes^r— . The data of the numerical (computer sim¬ 
ulation) and in situ experiments may be compared and 
understood in the framework of unified viewpoint only 
taking into account the realistic potentials of the inter¬ 
atomic interaction^. In addition to the physical nonlin¬ 
earity (which reflects the nonlinear character of the in¬ 
teratomic forces), the nonlinearity of the geometric ori¬ 
gin should be allowed for the framework of commonly 
used continuum approach to the CNT mechanics (the 
elastic thin shell theory, the beam theory or the finite 
element method)^. Besides the obvious consequncies 
of nonlinearity, such as a manifestation of the nonlin¬ 
ear elasticity or the loss of the stability under the large 
deformationi^ - — , there are some nonlinear effects which 
associate the CNT with the wide class of the nonlinear 
one-dimensional lattices^. It means that the effects of 
anomalous diffusion and the heat conductivity may be oc- 
cured in the CNTs^. The phenomenon of the slow energy 
exchange and the capture of the vibrations in some do¬ 
main of the CNT may appear in the CNT vibration s 17 ! 18 . 

The radial breathing mode (RBM) belongs to the most 
well known Raman-active oscillation branch of the car¬ 
bon nanotubes (CNT) vibrations and there are many 
studies of the RBM, which are based on the differ¬ 
ent approaches - since the continuum shell theory up 
to the quantum ab initio methodsiSr— 2 . The measure¬ 
ment of the RBM frequency allows to identify the CNT 
raduis22r— as well as to determine the critical pres¬ 
sure of the structural transitio n 26 ! 27 . The role of radial 
breathing vibrations for the transport properties of wa¬ 
ter molecules across SWCNTs has been also discussed 
recently^. One should note that the most of the CNT 
dynamics studies deal with the molecular dynamics (MD) 
simulation or the continuum aproach (the elastic thin 
shell or beam theory). However, because of the extreme 


complexity of the phonon band structure the dynamics 
of the different types of the CNT vibartions can not be 
described in the framework of the single approach. In 
particular, the reduced nonlinear theory of the elastic 
thin shell under appropriate physical hypothesis is re¬ 
quired to study the low-frequency circumferential flexure 
mode s 18 ! 29 . To investigate analytically the nonlinear ef¬ 
fects, which may occur in the radial breathig branch of 
the CNT vibration, a new model in the framework of the 
thin elastic shell theory is needed. 

However, one should remember that the RBM belongs 
to the optical branch and the specific approach is re¬ 
quired to analyse the nonlinear effects in the long wave 
domain, in view of the eigenfrequencies crowding and, 
as a sequency, a possibility of the resonant interaction 
of the nonlinear normal modes (NNMs). The adequate 
study of such highly non-stationary resonance processes 
is impossible in the framework of the NNM concept be¬ 
cause of the closeness of eigen frequencies that leads to 
appearance of a slow time scale and to the strong mutual 
interaction of the resonating NNMs. The key event in 
the stationary dynamics of the systems like the Fermi- 
Pasta-Ulam, sine-Gordon lattices is instability of one of 
resonating NNMs, which can be treated as an origin of 
the vibration localizatio n 30 ! 31 . It was shown recently^ 
that there is a need an alternative approach to describe 
the highly non-stationary resonance processes. This ap¬ 
proach is based on the new dynamical notion which is 
connected with the concept of Limiting Phase Trajec¬ 
tory (LPT) describing the maximum possible energy ex¬ 
change between different domains of the nonlinear sys¬ 
tem. The NNMs can be considered as the ground states 
in the slow time, and the motions, which are generated 
by their resonance interaction are new elementary exci¬ 
tations. Such consideration is the convinient when the 
initial conditions provide the excitation only the resrict 
set of the NNMs with close frequencies (the low temper- 
atutes and/or the selective initiations of the system os¬ 
cillations). The most intuitive example of such system is 
a pair of weakly coupled nonlinear identical oscillators^. 
Instead of conventional transition to modal presentation 
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(in the terms of in-phase and out-of-phase NNMs) the de¬ 
scription in the terms of the oscillators themselves turns 
out to be more adequate, because on the contrary to the 
resonating modes they are weakly interacting. As it was 
shown early, similar behaviour is the specific one in the 
quasi-one-dimensional nonlinear extensive systems (non¬ 
linear lattice s 17 ! 33 ! 34 and CNTs^). Some domains of the 
mentioned systems may be singled out, the behaviour of 
which turns out to be similar the nonlinear oscillators 
and the energy may be exchanged between them or may 
be localized in one of them. 

In this paper we consider the resonant interaction of 
the nonlinear normal modes (NNMs) near the left (low- 
frequency) edge of the radial breathing branch in the 
framework of the nonlinear dynamical equation for the 
radial component of the displacement field. This equa¬ 
tion was obtained from the nonlinear Sanders-Koiter thin 
shell theory in the asymptotic limit under hypotesis of 
the smalness of the Poisson ratio and it likes the non¬ 
linear Schrodenger equation with the specific type of 
the nonlinearity. The spectra of the RBM branch for 
the CNTs with various aspect ratios and under different 
bondary conditions are compared with those obtained by 
the direct numerical integration of the equations of the 
Sanders-Koiter thin shell theory. The dynamical regime 
of the intensive energy exchange results from the nonlin¬ 
ear normal modes interaction. The direct numerical inte¬ 
gration of the nonlinear equations of the Sanders-Koiter 
theory confirm the estimations following from the ana¬ 
lytical model. 



FIG. 1. The CNT spectrum according to the exact Sanders- 
Koiter thin shell theory: solid curves correspond to circum¬ 
ferential wave number n — 0, dashed ones - to n = 1 and 
dot-dashed one - to n = 2. The insert shows the small wave 
number part of the CFM branch. All the frequencies u are 
measured in dimensionless units and k - denotes the number 
of longitudinal half-waves along the CNT. 


II. THE MODEL 


The long wave length dynamics of carbon nanotubes 
is the obvious area of applicability of the classical the¬ 
ory of thin elastic shell, because the restrictions resulting 
from the plastic deformation are absent at the nanoscale 
level. The only complicating factor is the uncertainty of 
the parameter characterizing the thickness of the CNT—. 
However, this parameter does not play any significant 
role in the description of the axisymmetric deformations 
like the radial breathing vibrations. In such a case the 
momentless thin shell theory is a well approximation for 
the description of the long wave length RBMs. As it will 
be shown further, even the nonlinear effects don’t change 
this conclusion. The applicability of a well-designed thin 
elastic shell theory allows us to obtain an effective de¬ 
scription of the vibrational spectrum in the framework 
of the linear approximation. 

It is convenient to use the dimensionless variables 
which determine the elastic deformation of thin circular 
shell. In such a case all components of the displacement 
field (u - longitudinal, along the CNT axis, v - tangential 
and w - radial, respectively) are measured in the units 
of the CNT radius R. The displacements and respective 
deformations refer to the middle surface of the shell. The 
coordinate along the CNT axis £ = x/L is measured via 


the length of nanotube and varies from 0 up to 1, and 
ip is the azimuthal angle. One can define the dimension¬ 
less energy and time variables, which are measured in the 
units Eq = YRLh /(1 — v 2 ) and to = 1/ y/Y/pR 2 (l — 
respectively. Here Y is the Young modulus of graphene 
sheet, p - its mass density, v - the Poisson ratio of CNT, 
and h is the effective thickness of CNT wall. There are 
two dimensionless geometric parameters which character¬ 
ize CNT: the first of them is inverse aspect ratio a = R/L 
and the second - effective thickness shell (3 = h/R. 

The energy of elastic deformation of the CNT in the 
dimensionless units is written as follows: 
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where e^ and are the longitudinal, circumferen¬ 
tial and shear deformations, and and are the 

longitudinal and circumferential curvatures, and torsion, 
respectively: 
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Taking into account that the radial breathing normal 
modes do not depend on the azimuthal angle ip (the re¬ 
spective azimuthal wave number n = 0) and the transver¬ 
sal displacement v == 0, one can write the equation of 
motion as follows: 



FIG. 2. (Color online) The RBM spectra for the CNTs with 
different aspect ratios according to eq. © (solid lines) and 
eq. @ (dashed lines). The aspect ratios are equal to 20 (blue 
and dashed yellow), 30 (green and dashed orange), and 40 
(violet and light braun). 


Taking into account the expression © one can find the 
relation between u and w in the coordinate space as 


d 2 u 


dw d 2 w 
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Dispersive curves of the linearized problem contains two 
branches: 


w 2 = ^ ^1 + a 2 k 2 db y^l — ak ) 2 + 4<a 2 z ; 2 k 2 ^ (5) 
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Using this ratio one can rewrite second equation of © 
as follows: 
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To perform the asymptotic analysis of the long wave¬ 
length dynamics of the RB modes, it is convenient to 
rewrite equation © using complex variables: 


the one of which corresponds to the longitudinal acous¬ 
tic modes (sign ”+”) and the second one - to the radial 
breathing modes (sign ”-”). k is the longitudinal wave 
number. Figure [2] shows the eigenfrequencies of the long 
wavelength modes for the CNTs with different aspect ra¬ 
tios. In spite of that the real wave number for the peri¬ 
odic boundary conditions starts from the value k = 1, the 
product ak is small if the CNT is long enough. The long 
wavelength limit of © shows the spectrun crowding: 
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Returning to expression © one can see that the as¬ 
sumption of the Poisson ratio smallness v < 1 allows to 
use it as a small parameter. Using the multiple time 
scales ro = r, r\ — i zr, T 2 = u 2 t, etc., and expanding the 
function ip into series of the small parameter v\ 


w ~ 1 + i a 2 u 2 k 2 (6) 

The respective eigenvector 

(u, w) — (—auk, 1) (7) 

shows the relationship between longitudinal and radial 
components of the displacement field. 


ip = V (ipo + uip\ 4- U 2 ip 2 + ... ) e ZT , 

one can get the equation for the main order amplitude 
in the ’’slow” time i ~2 (see, for example j 17 i 36 for details): 
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Equation m admits the plane-wave solution 


Ipo = A exp (-i (wr 2 - &£))• (11) 

Taking into account the ” slowness” of the time 72 , the 
respective dispersion ratio should be written as follows: 

UJ = 1 + V 2 UJ = 1 + a 2 v 2 k 2 + -a A u 2 k A A 2 . (12) 
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One can see that the amplitude-independent part of 
m is in accordance with the relation (J6j), while the ef¬ 
fective positiveness of nonlinear addition points out the 
hard type of nonlinearity. 

Equation (HOl) is the modified Nonlinear Schrodinger 
Equation (NLSE) with the gradient type of the nonlin¬ 
earity. As it is well known, the standard NLSE admits the 
localized solution - the envelope soliton or the breather. 
However, a possibility any localized solutions in equation 
m is unknown. We will try to examine the possibility 
of energy localization while dealing equation m- First 
of all we replace equation (HOl) by its modal representa¬ 
tion, taking into accout only two resonant NNMs with 
the wave numbers k\ and 


V’o = Xi(72)sin(7rfci£) + X 2 (r 2 ) sin (nk 2 £) (13) 


After substitution of solution m into equation m 
one should use the Galerkin procedure to obtain the 
equations for complex amplitudes xi and X 2 • 
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where dwi = ^7 r 2 & 2 , (i = 1,2) are the modal frequency 
shifts (in the ’’slow” time scale 72 ) from the boundary 
frequency cjq = 1 of the considered brunch and 


°i j = ^ a 4 k i k j (*>i = M). ( 15 ) 

One can see that the nonlinear terms in equations 
m are separated into two groups: the terms \Xj\ 2 Xi 
(i,j = 1,2) determine the nonlinear frequency shift, 
while the terms XiXj(i 7 ^ j) describe the nonlinear inter¬ 
action between modes. The Hamiltonian corresponding 
to equations m can be written as 


H = |xi | 2 +< 5 w 2 |X2 | 2 - ^ (o-ii |Xi | 4 + 022 |X2| 4 ) 

-0-12 ^2 Ixil 2 IX2I 2 + t (x?X2 2 + X^X 2 )) 

(16) 


Besides the obvious energy integral equations 

m possess another integral 


* = |xi | 2 + Ix 2 | 2 , (17) 


which characterises the excitation level of the system. 
This is an analogue of the occupation number integral in 
quantum-mechanical terminology. 

As it was noticed in 1 ^ 18 ! 36 ! 37 , the modal description 
of nonlinear dynamics turns out to be inconsistent under 
resonant conditions. Therefore we need in the introduc¬ 
tion of new weakly interacting variables We have shown 
earlyii that they are the linear combination of the res¬ 
onating modes: 


1 1 

0i = ^=(xi + X 2 ); 02 = —^=(xi~ X 2 )- (18) 


The new variables describe the dynamics of some parts 
of the CNT 18 (or some groups of the particles in the effec¬ 
tive discrete one-dimensional chai n 17 i 33 i 34 ), because one 
can see a predominant energy concentration in certain 
domain of the CNT, while the other part of CNT has a 
lower energy. Due to small difference between the modal 
frequencies, the selected parts of CNT demonstrate a co¬ 
herent behavior. They are similar to two weakly coupled 
nonlinear oscillators and their temporal behavior looks 
very like to the beating in such system. Therefore we 
can consider these regions as new large-scale elementary 
blocks, which can be identified as unique elements of the 
system - the ’’effective particles”—. 

One can notice that the introduction of effective par¬ 
ticles m makes equations to be more complicated. 
However, due to the presence of the integral the 

dimensionality of the phase space of the system can be 
reduced. The ’’occupation number” X parametrizes the 
total excitation of the system, but the distribution of 
the energy is determined by the amplitudes of the ’’ef¬ 
fective particles” as well as by the phase shift between 
them. Actually, taking into account expression m , one 
can describe the behavior of the ” effective particles” with 
two real functions: 


= VXcosOe zA / 2 ; 02 = Vx smOe lA / 2 , (19) 


where the variable 6 characterizes the relative ampli¬ 
tudes of the ’’effective particle” and the variable A - the 
phase shift between them. 

Substituting relationships mm into equations d, 
we obtain the equations of motion in the terms of ’’an¬ 
gular” variables (#, A): 
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FIG. 3. (Color online) Phase portraits of the system (l20l) 
for the aspect ratio L/R = 80 and the ” occupation number” 
X = 2.5. The LPTs mark out by the thick black lines. 
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There are two types of the fundamental solutions on 
the presented phase plane. The stationary points cor¬ 
responding (in the slow time) to NNMs determine the 
stationary dynamics of the system. However, other type 
of phase trajectories is significant for understanding and 
description of highly non-stationary resonant dynamics. 
These are the trajectories, which separate the NNMs at¬ 
traction domains and they are the most distant from the 
stationary points. Such trajectories correspond to the 
extremely non-uniform distributions of the energy (from 
a possible ones). They were classified as the Limiting 
Phase Trajectories (LPTs). The main distinction be¬ 
tween the LPT and the separatrix is that the period of 
the motion along them is, generally speaking the finite 
one. The motion along the LPT between the states 0 = 0 
and 0 = 7r/2 leads to the redistribution of the energy be¬ 
tween the different parts of the system, i.e., between the 
effective particles. This process is an analogue of the beat¬ 
ing in the system of weakly coupled oscillators. It was 
shown that an adequate temporal description of LPT can 
be obtained in the terms of non-smooth functions which 
are saw-tooth function and its derivative in the sense of 
the distributions theory. 

All of these peculiarities of the phase space of the sys¬ 
tem (j20j) are shown in Fig. [3] The representative do¬ 
mains of the phase space are bounded by the intervals 
0 < 0 < 7r/2 and — n/2 < A < 37r/2. 

The numerical solutions of Eq. ([20]) with the initial 
conditions corresponding to the immovable point (0 = 





FIG. 4. (Color online) Time evolution of the variables 6 (red 
solid lines) and A (blue dashed lines) for the CNTs with dif¬ 
ferent aspect ratios: (a) L/R — 20; (b) L/R —. 80. The 
”occupation number” X = 0.5. One should pay the attention 
that r is the ’’slow” time and the real times (in the own period 
of the RBM) result as the value of r divided by the square of 
the small parameter. 



FIG. 5. (Color online) Period of energy exchange vs inverse 
aspect ratio of the CNT. Blue, orange and green curves cor¬ 
respond to the excitation X = 0.01, 0.5, and 1, respectively 


7r/2, A = 7r/2) for two values of the excitation X are 
shown in the Fig. H](a-b). 

The solutions, which are shown in Fig. 0J correspond 
to a slow redistribution of the energy between the ” ef¬ 
fective particles”. If the initial conditions respect to the 
LPT, the energy exchange reachs the maximum of the 
possible amount. The period of such energy exchange 
may be estimated as the time of passing the trajectory: 

(2I » 

The integral in equation (l2T|) can be estimated from the 
first of equations (1201) taking into account that A ~ ±7r/2 
on the LPT (see Fig. [3]) and the transition from A = 
—7r/2 to A = 7 t/2 takes no time: 


T _ 2 r ' 2 de _ 

Jo d6/dT2 {{ 5 u >2 — £>Wi) + | (<X 22 — ^ll) X) 

( 22 ) 

The variation of the period ([22]) with the aspect ratio 
of the CNT is shown in Figj5] 

To analyse the possibility of the stationary states insta¬ 
bility one should rewrite Hamiltonian ([1 6]) in the terms 
of the variables 6 and A: 
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H( 6 , A) = — ((5ui + 6002 ) + (Swi — Suj 2 ) cos A sin 26) 

A 2 

+ -r7r(3<T22(l — cos Asin2#) 2 + 3o-n(l + cos Asin2#) 2 
16 

+ 4(7i2 (3 cos 2 26 + sin 2 2# sin 2 A)). (23) 

The conditions of instability may be formulated as: 
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at (6 = 7r/4, A = 0) or (0 = 7r/4, A = 7r). 
The latter results in 


(24) 


x= 2 A = Q 
3(<tii - 2cr 12 ) 
2 (^ 2 -^) 

A = ——- 7 A = 7T, 

3(2(j 12 — 0 - 22 ) 


(25) 


Taking into account the definition of Scjj and relation¬ 
ship m one can see that no bifurcation at the pos¬ 
itive values of ”occupation number” X exist. There¬ 
fore, no localized breather-like excitatons can exist in 
the single walled CNT. So, in the contrary with the low- 
frequency circumferential flexure vibrations^, only the 
intensive energy exchange is possible in the radial breath¬ 
ing branch. 

The analysis performed above is based on the asymp¬ 
totic expansion of the equations in the framework of the 
nonlinear Sanders-Koiter elastic thin shell theory. Only 
two modes in the RB branch were taking into account. 
Therefore, a verification of our conclusion by the inde¬ 
pendent numerical methods is needed. Fig. [2] shows that 
if the aspect ratio of the CNT is large enough, more than 
two modes can be under resonant conditons. So, the in¬ 
fluence of the other part of the spectrum is very impor¬ 
tant for the estimation of the reliability of the obtained 
results. Our approach consists in the direct numeri¬ 
cal integration of the modal nonlinear equations of the 
Sanders-Koiter thin shell theory. The detail procedure 
was described for the circumferential flexure modes in-22. 
Therefore, we consider the method extremely shortly. In 
order to carry out the numerical analysis of the CNT dy¬ 
namics, a two-step procedure was used: i) the displace¬ 
ment field was expanded by using a double mixed series, 
then the Rayleigh-Ritz method was applied to the lin¬ 
earized formulation of the problem, in order to obtain an 
approximation of the eigenfunctions; ii) the displacement 
fields are re-expanded by using the linear approximated 
eigenfunctions, the Lagrange equations were then con¬ 
sidered in conjunction with the nonlinear elastic strain 
energy to obtain a set of nonlinear ordinary differential 
equations of motion. 

To satisfy the boundary conditions the displacement 
field was expanded into series 


r(€,<p,t) = 


M u N 


]T]TiWC(0 


m=0 n=0 


m (26) 


where function r(£, ip, t) substitutes the displacements 
u , v or w. 

In equations ([26]) T^(£) = T m (2£ — 1) are the Cheby- 
shev orthogonal polynomials of the m — th order and f(t) 
describes the time evolution of the CNT vibrations. 

Such an expansion allows to estimate the natural fre¬ 
quencies (eigenvalues) and modes of vibrations (eigen¬ 
vectors) under various boundary conditions. The results 
of performed calculation show that the eigenspectrum 
values are in the good accordance with the estimations 
made in the framework of reduced Sanders-Koiter theory 
discussed above (see Fig. ED. 

In the nonlinear analysis, the full expression of the 
dimensionless potential energy E e i containing terms up 
to the fourth order (cubic nonlinearity), is considered. 
Using the Lagrange equations 


77+ 7^=°’ r = {u,v,w} (27) 

a set of nonlinear ordinary differential equations is ob¬ 
tained; these equations must be completed with suitable 
initial conditions on displacements and velocities. This 
system of nonlinear equations of motion was finally solved 
by using the implicit Runge-Kutta numerical methods 
with suitable accuracy, precision and number of steps. 
The solution of nonlinear equations with initial condi¬ 
tions in the vicinity of the the LPT shows the energy ex¬ 
change process, the period of which coincides with equa¬ 
tion ([22]) for the wide interval of aspect ratios and exci¬ 
tation amplitudes (see Fig. ED . 


III. CONCLUSION 

As a summary we would like to notice that the phe¬ 
nomenon of the partial or intensive full energy exchange 
is the specific one for the resonating nonlinear normal 
modes. The generalized form of the Hamiltonian (fl6) ) 
reflects the most common sutiation for nonlinear interac¬ 
tion in very different system a 17 d 8 i 32 i 33 i 37 . The result of 
this interaction is determined by the relations between 
the parametrs of nonlinearity (<Jij) and the frequency dif¬ 
ferences (ScOj). The most important notion in that is the 
Limiting Phase Trajectory, which corresponds to the ” el¬ 
ementary process” associated with the energy exchange 
or capture. One should emphasize that only the analysis 
performed in the framework of this conception is capable 
of predicting of the main bifurcation in the dynamics of 
the system, in spite of that the quantitative results may 
be deficient due to restricted description of the system. 

Concerning with the radial breathing vibrations the 
absence of energy capture and localization is rather un¬ 
expected result, because the instability of low-frequency 
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a) 


b) 


c ) 


FIG. 6. (Color online) The energy exchange in the CNTs with different aspect ratios: (a) L/R = 20, (b) L/R = 40, (c) 
L/R — 80. The initial ”occupation number” X — 0.5. The dark violet and light beige areas correspond to low and high density 
of the energy. 





breathing modes seems intuitively evident. However, the 
presence of the strong gradient nonlinearity in equation 
m breaks down this viewpoint. It is the bright distinc¬ 
tion of the radial breathing branch from the circumfer¬ 
ential flexure one^. 
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